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A brief outline of the Lorentz Integral Transform (LIT) method is given. The method is well 
established and allows to treat reactions into the many-body continuum with bound-state like tech¬ 
niques. The energy resolution that can be achieved is studied by means of a simple two-body 
reaction. From the discussion it will become clear that the LIT method is an approach with a 
controlled resolution and that there is no principle problem to even resolve narrow resonances in 
the many-body continuum. As an example the isoscalar monopole resonance of '*He is considered. 
The importance of the choice of a proper basis for the expansion of the LIT states is pointed out. 
Employing such a basis a width of 180(70) keV is found for the '*He isoscalar monopole resonance 
when using a simple central nucleon-nucleon potential model. 


I. INTRODUCTION 


In the last years ab initio approaches have gained an 
ever increasing importance in the calculation of bound- 
state and reaction observables of nuclear systems with 
a nucleon number A > 4 (for a comparison of various 
techniques see [3, i). A particularly challenging as¬ 
pect in such calculations is the proper treatment of reso¬ 
nances. Recently, via the EIHH expansion technique Q, 
the Lorentz integral transform (LIT) method [J, Q was 
applied to calculate the isoscalar monopole strength of 
the "^116(6, e') reaction with realistic nuclear forces pay¬ 
ing special attention to the role of the low-lying 0+ res¬ 
onance of "‘He While the resonance strength could 
be obtained, it was not possible to determine the width 
of the resonance. As it was explained the origin of this 
deficiency arises from the too low density of LIT states in 
the resonance region, but it was not really clear why the 
density happens to be so low. One may think of various 
possible reasons: (i) the LIT method is based on bound- 
state techniques and might be incapable to describe the 
shape of narrow resonances, (ii) the used expansion in 
hyperspherical harmonics (HH) was not extended to a 
sufficiently large basis, or (iii) the HH expansion itself 
is not able to have a sufficient density of LIT states in 
the resonance region. The first possibility can be safely 
ruled out. In fact in [3] it was shown that for a fictitious 
two-nucleon system with a narrow resonance in the ^Pi 
partial wave exact information about the shape of the 
resonance could be obtained with the LIT method. Thus 
the aim of the present paper is to study which of the 
two remaining reasons is re^onsible for the partial fail¬ 
ure of the LIT method in [6| and how one can improve 
the calculation for having a precise information on the 
resonance width. 


central NN potential model (TN potential) is used that 
had been previously employed in the very first LIT appli¬ 
cations for the calculation of the electrom^netic break¬ 
up of "‘He by electrons and real photons [1,Q. It is inter¬ 
esting to note that in @ one already finds a clear signal of 
the 0+ resonance in the "‘He longitudinal response func¬ 
tion, however, at that early stage of LIT applications the 
resonance was not studied in greater detail. 

The paper is organized as follows. Section H contains 
a brief description of the LIT method. In section HI the 
question of the density of the LIT states is addressed 
and illustrated for the above mentioned two-body case. 
A LIT calculation for the isoscalar monopole strength of 
the '‘He(e, e') is discussed in section IV, where it will be 
shown how an HH calculation can be modified in order 
to determine the width of narrow resonances. Finally, a 
summary is given in Section V. 


II. LORENTZ INTEGRAL TRANSFORM (LIT) 


Over the years the LIT approach has been applied 
to a variety of inelastic electroweak reactions. A rather 
large number of applications can be found in the review 
articles 0 , 0 . 

The LIT of a function R{E) is defined as follows 

L{a) = j dE C{E,a)R{E), (1) 


where the kernel £ is a Lorentzian, 


C{E,<7) 


I 

{E - anY + cr| 


( 2 ) 


To clarify the matter first a classical two-body prob¬ 
lem, namely deuteron photodisintegration in unretarded 
dipole approximation, is discussed. It is shown how the 
density of LIT states can be increased in this case. In a 
next step it is considered how one can achieve such an 
increase also for the above mentioned isoscalar monopole 
strength of the ^He(e,e') reaction. In order to do so a 


(tr = gr -\- icT/); the parameter <7/ controls the width of 
the Lorentzian. Because of the adjustable width, and dif¬ 
ferent from many other integral transforms, the LIT is a 
transform with a controlled resolution. However, aiming 
at a higher resolution by reducing cr/ might make it nec¬ 
essary to increase the precision of the calculation. This 
point will be discussed in greater detail in section HI. 
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For inclusive reactions the LIT L{a) is calculated by 
solving an equation of the form 

{H-a)^ = S, (3) 

where H is the Hamiltonian of the system under con¬ 
sideration and S is an asymptotically vanishing source 
term related to the operator inducing the specific reac¬ 
tion. The solution T is localized. This a very important 
property, since it allows to determine 'I' with bound-state 
methods, even for reactions where the many-body con¬ 
tinuum is involved. 

Having calculated d' one obtains the LIT from the fol¬ 
lowing expression 

L(ct) = (^|^). (4) 

The response function R{E) is determined from the cal¬ 
culated L(a) by inverting the transform. A general dis¬ 
cussion of the inversion and details about various inver¬ 
sion methods are given in iSini,!!!. 

An alternative way to express the LIT is given by 

L{a) = -^Im({S\ ^ 

(TJ V (J/j -j- 1(7J — 11 / 

This reformulation is useful since it allows a direct ap¬ 
plication of the Lanczos algorithm for the determination 
of L(cr) dll- In fact the calculations discussed in the fol¬ 
lowing sections are performed using expansions on basis 
sets with a subsequent use of the Lanczos technique. 

In order to calculate a specific reaction one has to spec¬ 
ify the source term S in Eqs. and ([S]) . In case of an in¬ 
clusive electroweak reaction response functions have the 
general form 

R{Ef) = J dEf\{fm\^SiEf -Eo-ui), ( 6 ) 


The LIT result then reads 


L{(j) = 


N 

E 

i— 


^ {an - EnY + a] ' 


(9) 


Note that this result is related to the so-called Lanczos 
response i?Lnczs by 

ALnczs(w, a-i) = —L{aR = Eq+lo, ct/) . (10) 

TT 

In the limit ct/ —>■ 0 the Lanczos response is equal to the 
true response function R. However, one often calculates 
Aliiczs for a small but finite ai value and identifies the 
Lanczos response with the true response, which in general 
is an uncontrolled approximation. In the LIT approach 
one does not make such an identification of transform 
and response function. A proper treatment requires an 
inversion (see discussion in Q). 

It is important to note that the definition of the LIT 
in Eq. (I) contains the full response function R with all 
break-up channels. For the calculation of the LIT one 
may use any complete localized A-body basis set. Au¬ 
tomatically for any such set strength from all break-up 
channels is contained in the LIT. However, in principle, 
it can happen that in a specific energy interval of a given 
reaction one basis set is more advantageous than another 
one. In fact, such a case is discussed in section IV. 


III. A SIMPLE TWO-BODY PROBLEM 

To better illustrate the LIT energy resolution a sim¬ 
ple two-body case is discussed, namely deuteron photo¬ 
disintegration in unretarded dipole approximation. The 
corresponding cross section is given by 


where |0) and |/) are ground-state and final-state wave 
functions of the system under consideration, while Eq and 
Ef are the corresponding eigenenergies and uj = Ef — 
Eq is the energy transferred to the system, finally, 9 is 
the specific transition operator that induces the reaction. 
Note that different from the normal convention here the 
argument of the response function R is Ef and not uj. 

The source term S for inclusive reactions has the fol¬ 
lowing form 

\S) = 0|O). (7) 

A solution of the LIT equation (O via an expansion 
on a basis with N basis functions can be understood as 
follows. One determines the spectrum of the Hamilto¬ 
nian for this basis thus finding N eigenstates (jjn with 
eigenenergies E„. The energies define the positions 
of the above mentioned LIT states. Furthermore, the 
LIT solution assigns to any eigenenergy a Lorentzian with 
strength S'„ and width cr/. It should be noticed that the 
source term \S) solely affects the strength leading to 

5n = |(^n|^|0)P- (8) 


<^unret(*^) — ^TT CX UJ Runreti^f — Rq ; (H) 

where uj denotes the photon energy and a is the fine 
structure constant. The transition operator for the re¬ 
sponse function i?unret Is the dipole operator, 

( 12 ) 

i=l 

where Zi and Ti^z are the z-components of the position 
vector and of the isospin operator of the ith nucleon, 
respectively. In case of the deuteron the dipole operator 
induces only transitions to the np final states ^Pq, ^Pi, 
and ^P 2 -^F 2 . For simplicity in the following example only 
transitions to the ^Pi partial wave are considered. The 
ansatz for the corresponding T reads 

|^)=i?(r) |(; = l,5 = I)j = I)|T=l), (13) 

where r, I, S', j, and T = 1 is the relative distance, orbital 
angular momentum, total spin, total angular momentum, 
and isospin of the np pair, respectively. The resulting 






FIG. 1. (Color online) Spectrum of the Hamiltonian eigenen- 
ergies for the ^Pi NN channel with the AV18 NN potential 
for the basis of Eqs. (I13I14II with various values of N and b: 
N = 10 with b = 1 fm, 0.5 fm, and 0.25 fm in (a) and b — 0.5 
fm with N = 10, 30, and 50 in (b). 

LIT equation can be easily solved by direct numerical 
methods or by expansions of R{r) on a complete set. For 
nuclei with A > 2 very often HH expansions are used with 
separate hyperspherical and hyperradial parts, where the 
latter is usually expanded in Laguerre polynomials 
times an exponential fall-off. Therefore a corresponding 
ansatz is made here for R(r), 

N 

i?(r) =r^c„L£i)(r/6) exp(-^) , (14) 

n—0 

where c„ are the expansion coefficients and 6 is a param¬ 
eter regulating the spatial extension of the basis. 

For the following LIT results of the ^Pi channel the 
AV18 NN potential 0 is used. First the Hamiltonian 
eigenvalues entering in Eq. ®, are studied for var¬ 
ious basis sets. In the upper panel of Fig. 1 results are 
shown with 11 basis functions {N = 10) and different 
values for the extension parameter b. One sees that a 
greater spatial extension of the basis functions leads to a 
shift of the spectrum to lower energies. As lowest (high¬ 
est) eigenenergies one finds 5.17 MeV (7424 MeV), 1.27 
MeV (2689 MeV), and 0.32 MeV (790 MeV) for b = 0.25, 
0.5, and 1 fm, respectively. To obtain a higher density 
of LIT states one has to increase the number of the basis 
states N. This is illustrated in the lower panel of Fig. 1, 
where the cases with N = 10, 30, and 50 are shown for 
b = 0.5 fm. It is evident that the increase of basis func¬ 
tions does not only lead to a higher density of LIT states, 
but also to an extension of the eigenvalues En both to 
lower and higher energies. In fact as lowest (highest) val¬ 
ues one has now 1.27 MeV (2689 MeV), 0.19 MeV (11056 
MeV), and 0.076 MeV (25457 MeV) for N = 10, 30, and 
50, respectively. If one chooses the energy range up to 
pion threshold one finds that in all three cases about 
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FIG. 2. (Color online) LIT L{a) of the ^Pi channel with 
ai = 10 MeV and b — 0.5 fm. (a): N — 10 and 50; (b): the 
ratio r of the LIT with N = 10, 20, 30, and 40 to the LIT 
with A = 50 (the values 0.995 and 1.005 are illustrated by 
full lines). 


two-thirds of the N LIT states are located therein. 

After having discussed the energy distribution of the 
LIT states we come to the transform L{a) itself. First, in 
Fig. 2, the case with ct/ = 10 MeV and b = 0.5 fm is illus¬ 
trated. The upper panel of the figure shows that an in¬ 
crease of the number of basis functions from 11 {N = 10) 
to 51 (A=50) does not lead to noticeable differences. In 
other words, in this case, the LIT is already quite well 
converged with a rather small basis. This nice agree¬ 
ment of both results is not obvious, since Fig. 1 exhibits 
rather different positions and densities of the LIT states 
for both cases. This has to be interpreted as follows. 
A discretization of the continuum has no direct physical 
meaning and leads in principle to random results. On the 
contrary one can use the discretization to calculate inte¬ 
gral transforms, as for example the LIT, which lead in 
convergence to a unique result despite of the randomness 
of the discrete eigenenergies. 

In order to find differences between the two results of 
Fig. 2a one has to present them in a different way, as 
is done in Fig. 2b. There one sees that the agreement 
is extremely good up to about 30 MeV and thus one 
can consider the LIT for ct/ = 10 MeV to be already 
converged in this energy range with only very few basis 
functions. Beyond 30 MeV the LIT is not yet converged 
with A = 10 and starts oscillating about the LIT with 
A = 50 with differences becoming even greater than 5%. 
Because of this regular oscillations about the more pre¬ 
cise result it should even not lead to serious inversion 
problems as long the response does not contain any spe¬ 
cific structure at higher energies. If one wants to check 
the existence of such structures one should improve the 
precision of the calculation by enlarging the basis. In 
fact with larger A values of 20, 30, and 40 one finds an 
increased high-precision range with relative differences 
compared to the case with A = 50 of less than 0.5% up 
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FIG. 3. (Color online) LIT L{a) of the ®Pi channel with 
ai = 0.1, 1, 2.5, and 10 MeV in all four panels (fe = 0.5 fm), 
but different N values in the various panels: 10 (a), 20 (b), 
30 (c), and 50 (d). 

to about 50, 60, and 85 MeV, respectively. However, in 
order to search for possible structures it is not the proper 
strategy to just enlarge the basis. One has to consider 
that a resonance with a width considerably smaller than 
CT/ is smoothed out in the LIT such that the details of 
the shape are hidden in tiny contributions to the trans¬ 
form. To disentangle the details one should reduce uj 
increasing in this way the energy resolution of the LIT. 

The situation for smaller values of aj is illustrated in 
Fig. 3, where LIT results are shown up to about pion 
threshold for cr/ = 0.1, 1, 2.5, and 10 MeV and for various 
values of N. For the smallest value, iV = 10 (Fig. 3a), 
one finds isolated Lorentzian peaks. In case of ct/ = 0.1 
MeV they appear already at low energy, for cr/ = 1 MeV 
at somewhat higher energy, and for tr/ = 2.5 MeV at even 
higher energy. Since the density of LIT states grows with 
growing N (see other panels of Fig. 3) isolated peaks are 
pushed to higher and higher energies if N is increased. 
In addition one notes that any decrease of cr/, i.e. any 
increase of the resolution, shrinks the convergence range 
for a given value of N. For the smallest ct/ value of 0.1 
MeV a rather strong oscillatory behavior is still present 
at very low energies even with iV = 50. In order to get 
a smooth and converged LIT also in this case one would 
need to increase N considerably. 

For a reliable inversion of L(cr) the transform has to 
be sufficiently converged for a given cr/. In particular iso¬ 
lated peaks of single Lorentzians should not appear, i.e. 
for any cr/j value one should have a significant contribu¬ 
tion from more than one Lorentzian. Of course, for the 
present two-body case one could increase N further with¬ 
out greater problems, but in a general many-body calcu¬ 
lation this might not be easily possible. On the other 
hand it is not necessary to work with a single cr/ value 
only. Considering again the LIT results of Fig. 3 with 
N = 50 one sees that for cr/ = I MeV a rather converged 
result is obtained up to about 10 MeV. Thus structures 


of a relatively small width would leave a visible signal 
in the transform in this energy range. In fact one can 
use this LIT in the low-energy range and combine it with 
a LIT with a larger cr/ for higher energies. In general 
one can proceed as follows, one defines a new transform 
T in an interval [cro,criv/] by combining M transforms 
Lm = L(cr//, cr/_m) such that only the LIT with the spe¬ 
cific cr/ = enters significantly in the energy range 

M-1 

F = Li/i-f- ^ /m—l))/mT^iMLM(l /m— l)- 

m=2 

(15) 

The function fm in Eq. m is a smooth cutoff from 1 
to 0 at cr/{ = aR^rn, for example fm{crR < <JR,m) = 1 
and fm{crR > cfR^m) = exp(-((cr/i - cr/{_^)/A)2"), where 
one has to take reasonable values for the parameters A 
and n. The coefficients am in Eq. (fT51) can be chosen such 
that size of the transform does not change too drastically 
from one energy range to the next. Thus, for the case of 
the LITs of Fig. 3 with N = 50, one could set cr/i = 1, 
cr /^2 = 2.5, and cr /_3 = 10 MeV with cr//^ = 10 and 
<^R .2 = 30 MeV. In order to improve the precision in 
the threshold region even further one could include the 
case with cr/ = 0.1 MeV, but one would need to further 
increase N. Alternatively, keeping N = 50, one could 
check the convergence behavior of a LIT with a somewhat 
larger cr/, for example cr/ = 0.25 MeV. 

The discussion above shows that the LIT is an ap¬ 
proach with a controlled resolution. In an actual calcu¬ 
lation one should check which is the lowest cr/ that leads 
in a specific energy range to a sufficiently converged and 
smooth LIT without that a single LIT state sticks out. 
Structures which are considerably smaller than such a 
a I value cannot be resolved by the inversion. A helpful 
criterion is given in Q (see discussion of Fig. 7 in a)- 

IV. ^HE ISOSCALAR MONOPOLE RESPONSE 
FUNCTION 

The isoscalar monopole response function M{q,Ef = 
Eq+uj) can be determined in inclusive electron scattering 
(g and w represent momentum and energy transferred 
by the virtual photon to the nucleus). In this case the 
transition operator 0 of Eq. ([6]) becomes g-dependent and 
takes the form 

^( 9 ) = joiqri ), (16) 

i=l 

where is the nucleon isoscalar electric form factor, 

Vi is the position of nucleon z, and jo is the spherical 
Bessel function of 0*^ order. 

Experimental investigations of the ^He(e, e') reaction 
in |15l - [l7l| revealed a O'*' resonance located less than 1 
MeV above the "^He break-up threshold with quite a nar¬ 
row width of about 250 keV. It is interesting to note 
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FIG. 4. LIT L{a) of M{q,uj) at q — 300 MeVVc (TN poten¬ 
tial) of '‘He with four-body CHH basis from [|. 


that very recently a EIHH-LIT calculation was carried 
out, where it is pointed out that the resonance might 
be interpreted as a breathing mode [l^. Unfortunately, 
this and the preceding EIHH-LIT calculations @ were 
only able to determine a resonance strength but not a 
resonance width since the density of LIT states in the 
resonance region was too low. On the other hand, with 
the experience made in the previous section it seems to 
be easy to increase the density of LIT states also in the 
region of the '^He isoscalar monopole resonance, namely 
by increasing the number of HH basis states. Unfortu¬ 
nately with an HH expansion this does not work very 
well in the energy region below the three-body break-up 
threshold. To illustrate this, the first LIT application Q, 
already mentioned in the introduction, is considered in 
the following. The calculation was carried out for inclu¬ 
sive electron scattering off '^He and used a correlated HH 
(CHH) basis, where NN short range correlations are in¬ 
troduced to accelerate the convergence of the HH expan¬ 
sion. The NN interaction (TN potential. Coulomb force 
included) consisted in a central spin-dependent interac¬ 
tion active only in even NN partial waves. In Fig. 4 we 
show an unpublished result for the LIT of the response 
function M from this calculation. One sees that there 
is only an isolated LIT state in the resonance region. 
In this specific case one finds the next LIT state only 
about 0.1 MeV below the four-body break-up threshold 
at cT/j = 0 MeV. On the contrary for positive values of an 
one has quite a high density of LIT states. An increase 
of the number of hyperspherical and/or hyperradial ba¬ 
sis states does not change the general picture of having 
only very few LIT states for negative an values^as in 
fact it was the case in the LIT-EIHH calculations [y, . 

In addition it was checked for the present work that an 
increase of the parameter b in the hyperradial functions 
(~ Ln\plb) exp(—/3/26)) in the four-body CHH calcu¬ 
lation of Q does not lead to any significant change con¬ 
cerning the low-energy density of LIT states. 

At first sight it is not understandable why an in¬ 
crease of the HH basis states does not have a signifi¬ 


cant effect on the density of LIT states in the energy 
range below the three-body break-up threshold. On the 
other hand one has to consider that the dynamical vari¬ 
able for a two-body break-up, i.e. the relative vector 
r 4 = r 4 — X]i=iU/3 between the free nucleon and the 
bound three-body system, does not appear explicitly in 
the HH formalism. Therefore it might be helpful to use a 
different basis, where is taken into account as variable 
explicitly. In the following, the four-body calculation is 
switched to such a basis. The new basis states consist 
in a product of a three-body CHH basis state for the 
first three nucleons times a single-particle basis state for 
the fourth nucleon. The CHH basis state is then again 
a product of a hyperspherical basis state V[ifa] (U) times 
a hyperradial basis state Rn{p), where the former de¬ 
pends on the grand-angle Q and is characterized by the 
grand-angular quantum number K and a set of other 
quantum numbers [a], whereas the latter is proportional 
to Ln\p/b) exp{—p/2b) (for a detailed definition of the 
HH basis see for example i). 

The single-particle states are defined as follows 

l®(i’4)) = X! \krnii) \sims^) , (17) 

where k, S 4 = 1/2 and = 1/2 are the orbital angular 
momentum, spin, and isospin quantum numbers, respec¬ 
tively, and TO/ 4 , rus^ and nit^ denote the corresponding 
projections. The radial wave function (/liki) is given by 

3 A^4 —1 ^ f 

= (?’!)'" n ■^(''*4) Cn^L^//}{r'Jbi) exp ( ^ j 

i^l 714^0 ^ ' 

(18) 

where /(r/ 4 ) = /(|r/ — r 4 |) is the NN correlation func¬ 
tion (same correlation as in CHH basis). The CHH basis 
provides antisymmetric states for the first three nucle¬ 
ons, whereas the product of the CHH basis states and 
the single-particle states has to be antisymmetrized in 
order to have totally antisymmetric basis states for all 
four particles. In order to calculate the '‘He ground state 
or the transitions induced by the action of the operator 
0(g) of Eq. (ITCl) on the ^He ground state one has to cou¬ 
ple the CHH basis states with those of Eq. dia to a total 
angular momentum L equal to zero, and also the total 
spin (isospin) wave function of the four nucleons has to 
be coupled to a total spin S (isospin T) equal to zero. 

The calculations described in the following proceed in 
the same way as it was carried out in [^, |^ , namely by 
using a nine-dimensional Monte Carlo integration for the 
evaluation of the various Hamiltonian and norm matrix 
elements. Such an approach leads to reliable results (see 
benchmark test 0 ). 

In Table U the convergence of the '‘He binding energy 
is illustrated. It is evident that the dominant contribu¬ 
tion comes from the basis functions where the three-body 
CHH basis is in a symmetric state and the relative angu¬ 
lar momentum U of the single-particle motion is equal to 
zero. Mixed symmetric CHH states together with ^4 = 0 



TABLE I. Convergence of ^He binding energy BE with the 
TN potential (Coulomb force included): U denotes the or¬ 
bital angular momentum of the single-particle motion (see 
Eqs. (11711811 1. K 3 is the grand-angular quantum number of the 
three-body CHH states, “sym” and “mixed” indicate that the 
CHH state is symmetric and mixed symmetric, respectively 
(note that for Z 4 > 0 two symmetric CHH states and only the 
lowest mixed symmetric state are taken). 


u 

K 3 

CHH symmetry 

# CHH states 

BE [MeV] 

0 

0 

sym 

1 

28.66 

0 

4 

sym 

1 

30.11 

0 

6 

sym 

1 

30.56 

0 

8 

sym 

1 

30.67 

0 

10 

sym 

1 

30.77 

0 

12 

sym 

2 

30.82 

0 

14 

sym 

1 

30.85 

0 

2 

mixed 

1 

31.28 

0 

4 

mixed 

1 

31.32 

1 

1,3,5 

sym and mixed 

3 

31.39 

2 

2,4 

sym and mixed 

3 

31.41 


enhance the binding energy by about 0.5 MeV. The con¬ 
tribution of states with ^4 > 0 is very small and leads to a 
further increase of about 0.1 MeV (antisymmetric CHH 
states have been neglected). The final result of 31.41(5) 
MeV agrees very well with the converged result of the 
four-body CHH calculation of 31.40(5) MeV. 

For the calculation of the LIT of the ^He isoscalar 
monopole response function M(q,Lu) hyperradial and ra¬ 
dial basis states are chosen with a rather large extension 
in space. This should lead to a sufficiently high density 
of LIT states at low energy (see Fig. 1). The following 
choice is made: 6 = 4 fm (hyperradial CHH states) and 
&4 = 1.33 fm (radial single-particle states). For the hy¬ 
perradial part 15 basis functions are taken, whereas the 
number of the radial single-particle basis states is kept 
variable and denoted by N 4 . Since the main interest of 
this investigation is concentrated on the low-energy part 
of M{q, to) only those three-body CHH basis states which 
lead to a significant contribution to the ^He binding en¬ 
ergy (see Table |T| are taken into account: symmetric 
CHH states up to K 3 = 6 and the mixed symmetric CHH 
state with K 3 = 2. All the remaining states contribute 
only with 0.42 MeV to the "^He binding energy. Thus it 
is reasonable to expect also a shift of the position of the 
“^He O’*' resonance by about 0.5 MeV towards higher en¬ 
ergies with respect to the four-body CHH calculation of 
Fig. 4 leading to a value of about -5.9 MeV. To take into 
account only four hyperspherical three-body CHH states 
reduces the numerical effort of the calculation quite a 
bit, but even with such a reduced basis the calculation 
requests a considerable numerical effort (note that due 
to the rather weak fall-off of the hyperradial/radial basis 



FIG. 5. (Color online) Spectrum of the Hamiltonian eigenen- 
ergies for the "^He(J’^ =0'^) states with the TN NN potential 
for the basis described in the text with A 4 basis functions for 
the radial single-particle basis. 

functions the relevant integration volume becomes much 
larger than for the bound-state calculation). 

In Fig. 5 the energy distribution of LIT states is shown 
up to the four-body break-up threshold with three differ¬ 
ent N 4 values. As it has been anticipated the density 
of LIT states increases with a growing number of radial 
single-particle states. Due to the interplay of CHH basis 
states and radial single-particle basis states the pattern 
is not as regular as in the two-body example of Fig 1. 

It is obvious that a further increase of N 4 would lead to 
an even higher density of low-energy LIT states, but, un¬ 
fortunately, the precision of the nine-dimensional Monte 
Carlo integration is not sufficiently high and a solution 
with N 4 > 20 becomes problematic. In this respect it is a 
drawback to work with a correlated HH basis because one 
looses the orthogonality. Nonetheless, as shown in the 
following, the present calculation allows a determination 
of the width of the ^He isoscalar monopole resonance. 

In Fig. 6 the LIT of M(q,u!) is shown for N 4 = 20 
and N 4 = 22. In the latter case the basis function with 
U 4 = 20 (see Eq. (fTBll ) is dropped because it makes the 
numerical solution of the corresponding eigenvalue prob¬ 
lem problematic (precision of Hamiltonian and norm ma¬ 
trices is not sufficiently high). The figure illustrates that 
with increasing resolution, i.e. with decreasing < 7 / the res¬ 
onance becomes more and more pronounced against the 
background. One also hnds that the anticipated peak 
position of about -5.9 MeV is roughly confirmed. 

In order to determine a resonance width the LIT has to 
be inverted. The procedure of the inversion in presence of 
a resonance is described in Q . As rule of thumb one can 
say that the chosen aj should not be much larger than 
the resonance width. However, in principle, if the LIT 
is calculated with a very high precision one can choose 
also a considerably larger aj. In fact in the model study 
0 it is nicely demonstrated that the obtained width is 
independent from the used <ti over a very large range, 
even with ai = 5 MeV the shape of a resonance with a 
width of 270 keV was exactly determined. In the present 
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FIG. 6. (Color online) LIT L{a) oi M{q,uj) at g = 250 MeV/c 
(TN potential) of "^He with ai = 0.01, 0.1, 0.25, 1, and 5 MeV 
calculated with the new basis states (three-body CHH basis 
times single-particle basis of Eqs. (I17I18II 1: in (a) the single¬ 
particle basis with 20 and in (b) with 21 states (see text). 


case the precision of the calculated LIT is not as high as 
in 0. The actual inversion was made with two different 
aj values (see discussion of transform T introduced in 
Eq. (HH)), (T /^2 = 5 MeV beyond -4 MeV, whereas 
was varied in the range from 0.1 to 0.5 MeV in the low- 
energy region. The various results for the width were 
quite stable (maximal difference: 0.01 MeV) and lead to 
the following values: 120(10) keV {N^ = 20) and 240 
keV {N 4 = 22), the mean value amounts to 180(70) keV. 
The result lies in the same ballpark as the experimental 
value of 270(50) keV [l^. One could try to increase the 
theoretical precision of the determination of the width 
by a further increase of N 4 . Because of the problem just 
mentioned above this would require a non negligible effort 
in the present calculation. On the other hand it is much 
more desirable to make such a calculation with a realistic 
nuclear force instead of having a very precise result for a 
simple NN interaction like the TN potential model. 

For the inversion of the LIT it is assumed that there is 
only a single peak in the resonance region. In principle 
from an increase of the precision of the calculation one 
could also find out that the resonance has a more compli¬ 
cated structure, for example a double peak. It is evident 
from Fig. 6 that in the present calculation one controls 
the resonance with a resolution of about 100 keV, thus it 
is not possible to resolve structures with an even smaller 
width. The situation is similar to experiment, where one 
works with a given resolution of the experimental appa¬ 
ratus. 


V. SUMMARY 

In this work first a brief outline of the LIT approach is 
given, a method that allows to calculate continuum ob¬ 
servables with bound-state techniques. In a next step the 


energy resolution that can be obtained with the method 
is discussed. To this end a specific channel of the final 
state in deuteron photodisintegration is studied in unre¬ 
tarded dipole approximation. In order to do so the LIT 
equation is solved using a bound-state basis to calculate 
the relevant quantities. It is shown how the spatial ex¬ 
tension and the number of basis states affects the calcu¬ 
lation. The role of the width parameter aj of the LIT is 
discussed, especially which aj values should be chosen in 
a specific situation. The discussion makes clear that the 
LIT constitutes an approach with a controlled resolution. 
In particular it allows to determine the width of narrow 
resonances in the continuum. However, for LIT calcula¬ 
tion with an arbitrary many-body bound-state basis this 
is not always guaranteed. Usually narrow resonances of 
an A-body system in the continuum are located below 
the three-body break-up threshold. If one intends to 
study the resonance in a scattering state calculation it 
is mandatory to take into account the relevant dynami¬ 
cal variable, i.e. the relative vector between the two frag¬ 
ments. In spite of the LIT bound-state character this 
dynamical variable should also appear explicitly in a LIT 
calculation. It guarantees that the density of LIT states 
can be enhanced in the resonance region by increasing 
the number of those basis states that directly depend on 
the dynamical variable. On the contrary, if this variable 
is not included explicitly, as for example in an A-body 
HH calculation, it is difficult, maybe even impossible, to 
obtain a detailed information about the resonance. 

In a general case of a resonance with various open chan¬ 
nels the resonance can have a partial decay width to all 
these channels, which then results in a total decay width. 
In this case it is sufficient to choose just one of the various 
possible dynamical variables. The only aim which has to 
be fulfilled is a sufficient density of LIT states. The LIT 
then by definition collects strength from all open chan¬ 
nels and a determination of the width via inversion leads 
to the total width. 

In order to illustrate the situation in greater detail the 
isoscalar monopole response function M of ^He is consid¬ 
ered as test case using a simple spin-dependent central 
NN interaction (Coulomb force is included). In fact, here 
one finds a narrow resonance in the continuum below the 
three-body break-up threshold. Obviously, in this case 
the dynamical variable is given by the relative vector of 
the free nucleon and the bound three-nucleon system. Us¬ 
ing a proper basis as described above it is indeed found 
that the density of LIT states grows in the resonance 
region if the number of single-particle basis states is en¬ 
hanced. It is shown that the LIT state density becomes 
sufficiently high to determinate the resonance width via 
an inversion of the transform. In fact, a width of 180(70) 
keV is found, a result that is not too far from the experi¬ 
mental value of 270(50) keV. It would be very interesting 
to perform such a calculation also with modern realistic 
nuclear forces as have been used in the LIT-EIHH cal¬ 
culations of the response function M[^, [l^, however the 
effective interaction approach has to be a bit redesigned, 





















since one would not have any more a pure HH basis. 

The present approach is not only advantageous in case 
of cross sections with narrow resonances, but also for non¬ 
resonant two-body break-up cross sections at very low 


energies. The possibility to increase the LIT state density 
allows to work with smaller aj values thus enhancing 
the energy resolution of the calculation. This might be 
particularly interesting in case of astrophysical reactions. 
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